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SUMMARY 


In an attempt to better predict turbulent-boundary-layer develop- 
ment in aircraft diffusers, a procedure based on the kinetic-energy 
equation and an extended form of the momentum equation has been devised 
for calculating the development of turbulent boundary layers, in adverse 
pressure gradients, for that class of flows for which the fluid density 
at all points and total pressure outside the boundary layer are 
invariant. In the development of this procedure an effort was made not 
only to arrive at an analytical form that would allow examination of 
the significance of the physical quantities involved but also to achieve 
a high degree of consistency with the more recent results of turbulent- 
boundary-layer research. 

Correlations, which are essential to the execution of the method, 
of several elements of the basic equations in terms of quantities 
regarded as controlling parameters are presented. These correlations 
are strictly interim empirical relationships both as to numerical values 
and as to the variables involved. It is to be expected that, with the 
acquisition of more data covering a broader range of conditions, improve- 
ments will be made in both the functional nature and the accuracy of the 
correlations. A comparable and associated refinement of the basic 
equations is likewise anticipated as more data are obtained. 

Predictions, by this method, of turbulent-boundary-layer develop- 
ment are compared with experimental results from several sources for a 
number of cases of flow over flat plates and airfoils and in conical 
diffusers. In the range of boundary-layer flow short of separation, the 
agreement of calculated values with experimental values was, in most 
cases, quite satisfactory; in some instances, however, definite disagree- 
ment between the calculated and experimental resxilts was noted. It is 
believed, however, that good agreement has been obtained in enough 
instances to Justify continuation of effort along the present lines, 
particularly with respect to improvement of the correlations and refine- 
ment of the equations. 
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INTRODUCTION 


Because of the dominant part played ly the boundary-layer growth 
in determining diffuser performance, theoretical diffuser analysis is 
effectively a problem in the calculation of the boundary-layer develop- 
ment’. For aircraft application the boundary layer is usually turbulent 
and of a thickness less than the half-span of the channel. A semi- 
empirical procedure for the calculation of the development of such 
boundary layers in two-dimensional flows has been given by Von Doenhoff 
and Tetervin (reference 1), Attempts at direct application of this 
and other existing methods to conical-diffuser flows, however, have led 
to difficulties both in the accuracy obtained under some conditions and 
in the physical interpretation of ttie basic equations. The present 
research was undertaken with the object of evolving a procedure which 
would facilitate inteipretation of the physical quantities and would be 
applicable to conical-diffuser flows. This present paper is primarily 
concerned with the. progress to date in achieving this objective. It 
deals with a procedure developed for that class of flows for which the 
fluid density at all points and total pressure outside-the boundary layer 
are invariant, presents interim en?)irical correlations of certain 
physical quantities essential- to the execution of the procedure, and 
includes an examination of the results obtained to date. 


SYMBOLS 


D duct diameter 

H boundary-layer shape parameter (6*/0) 


q 


dynamic pressure outside boundary layer 



Rq boundary-layer Reynolds number based on momentum 

thickness (9Up/[Ji) 

u velocity parallel to surface at a perpendicular distance 

y from wall 

u' root-mean-square fluctuating velocity in axial direction 

U velocity parallel to surface at a distance greater than 

or equal to 6 from wall 

X longitudinal distance parallel to wall measured from 

initial station 



NACA TN 2 U 78 


3 


^ X 


duct-station locations measured longitudinally on duct 
axis from initial station 


T 






c 


airfoil chord 




mean inlet Mach number 


7 

5 

5 * 

0 

P- 

P 

T:o/2q 

(Tc/2q)g^Y. 

(V2q)L.T. 

'^R/2q 

''^s/2q 

Subscripts : 
exp 


perpendictolar distance from wall to point at which 
velocity u is measured 

boundary-layer thickness, defined as perpendicular 
distance from wall to point at which contribution 
to integrals for 6* and 0 is negligible 



viscosity 

density 

wall shear-stress coefficient 


wall shear-stress coefficient obtained by equation of 
Squire and Young (reference 2) 

wall shear-stress coefficient obtained by equation of 
Ludwieg and Tillmann (reference 3) 


turbulent-normal-stress coefficient 


dissipation coefficient 



.: l . _2. 

2q 3y 



from experimental data 




corr 


from correlation curve 
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ANALYSIS 

Momentum Thickness 


In the two-dimensional analysis- of reference 1, the rate of growth 
of the momentum thickness was calculated by the momentum equation 


M = ZO _ i da/i + H\ 

dx 2q q dx\ 2/ 


( 1 ) 


The Squire and Young (reference 2) values of., skin-friction coefficients 
were used with this equation to calculate the rate of growth of the 
momentum thickness. 

The sin^slest apparent procedure for extending equation (1) for 
application to conical diffusers is to redevelop it in circular measure 
by simple summation of forces in the longitudinal direction as given 
in reference U, for example. If the assumption is made that the 
boundary-layer thickness is small in relation to the duct diameter, 
the following approximate equation results: 


^ _ i i ^ 

dx 2q q dx\ 2/ D dx 


This equation differs from equation (1) by the presence of an additional 
term which provides for the effect of the changing duct perimeter. The 
same equation can be derived more rigorously by starting directly with 
the Navier-Stokes equations. 

In general, the values of momentum thickness calculated 
equation (1) for two-dimensional flow and by equation (2) for conical- 
diffuser flow closely approximate the experimental values. Discrep- 
ancies, local in character, between the calculated and experimental 
gradients were, however, noted. Figure 1 shows typical examples of 
experiittental and calculated momentum-thickness gradient and momentum 
thickness as a fxmction of x as illustrated for the two-dimensional 
flow data from reference This figure shows that over most of the 
range quite good agreement is found between experimental and calculated 
momentum- thickness gradient and momentum thickness. Near the downstream 
end of the surface, however, the difference between experimental and 
calculated momentum-thickness gradient is appreciable; this difference 
is responsible for the gradual divergence of the two curves for momentum 
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thickness. Several investigators (references 6 to 9) have observed^ 
and commented upon discrepancies of this type, the occurrence of which 
seemed to be associated with proximity to separation. Newman (refer- 
ence 10) suggested that this lack of agreement was due to evaluation of 
the momentum equation in terms of mean velocity and the neglect of terms 
arising from the presence of velocity fluctuations. Wallis (refer- 
ence 11) mentioned that the term arising from the Reynolds normal stress 
is of increasing magnitude as separation is approached and implied that 
this term might eliminate these discrepancies if included in the momentum 
equation; inclusion of this term gives 


M = Is + _ i Sa/i + s') (3) 

dx 2q 2q q dxV 2/ 

in which the turbulent-normal-stress coefficient Tj^2q is given by 
the equation 




l£ 

2q 


r® 1 a(pu'2) 

0 2q 3x 


dy 


(h) 


The corresponding momentum equation for conical diffusers is 


^ = 12 + IS _ i 2a/i + S') _ i 22 

dx ~ 2q 2q q dx\ 2/ D dx 


( 5 ) 


An examination of equation (3) was made ty using the two-dimensional 
data of reference 5 and the results are illustrated in figure 2. The 

solid curve ('^ + — | was obtained ly solution of equation (3) in 
\2q 2qJ 

J A 

which the slopes ^ used were obtained from a curve faired through 

the experimentally measured momentum thicknesses. A curve of the Squire 
and Young skin-friction coefficient (To/2q).s,Y. given in this 

figure. The decreasing value of the Squire and Young skin-friction 
coeff icimt" with increasing longitudinal distance is due to increasing 
boxmdary-layer Reynolds number. The values of turbulent-normal-stress 
coefficient evaluated from hot-wire measurements in refer- 

ence 5 are included in figure 2. It is apparent that the sum of the 
turbulent-normal-stress coefficient r^.2q and the skin-friction 
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coefficient ^To/2q^g Y exceeds the value of the sum obtained by 

evaluation of equation (3) from the l8-foot station to approximately 
the 22-foot station. 

Experiments by Ludwieg and Tillmann (reference 3) Tqy a hot-plate 
technique have indicated that the skin-friction coefficient may diminish 
with increasing shape factor Hj a variable consideration of vAiich was 
beyond the scope- of the Squire and Young representation (reference 2). 
Values of skin-friction coefficient (To/2q)j^ ^ calculated from the 

Ludwieg and Tillmann correlation (reference 3) are also shown in figure 2. 
The sum of the values of skin-friction coefficient (xQ/Zq)-^ Y 

the turbulent-normal- stress coefficient Tj^2q gives a good approxi- 

Tq TJJ 

mation to the sum ^ ^ evaluated from equation (3) from the 

18-foot to the 22-foot station. The approximation ceases to be good, 
however, beyond the 22-foot station, the points departing greatly but 
without apparent system from the curve. This deviation may indicate 
significant incompleteness of equation (3) in this region. However, it 
may be noted that the ^Tp/2q)g^ values also depart greatly and 

apparently without system from the trend that might have been expected 
on the basis of preceding values. This departure suggests the alternate 
possibility that the origin of the deviations may lie in increased 
difficulty in determining the experimental value of x^2q at the 
downstream end of the surface. 

The change in turbulent=normal-stress coefficient with increasing x 
is noted in figure 2 to be opposite in direction to that of the skin- 
friction coefficient ('’^o/2q]L^;ii^ over this range, and the sum of 

('t^o/2^)l.T. 3^8 T|^2q does not differ greatly from values of 

^To/2q^g^Y, the-first several feet; 

Because the pressure-gradient term in equation (3) usually makes a 
very substantial contribution to the rate of growth of momentum 
thickness, the differences resulting from the choice of eq^tion (1) in 
which ('’^o/2q)s^Y. used or equation (3) in which ('Ry 

'’^R/2q are used are not so great as might be surmised from examination 
of figure 2. A better perspective is possible by reference to fig- 
ure 1(a), in which the use of_equation (3) with (T^o/2q)L.T. TTi/2q 

yields a curve coincident throughout most of its length with the experi- 
mental curve shown by the solid line in comparison with the dashed line 
which was computed from equation (1) with ('To/2q)s ^Y. • The deviation of 
the curve obtained by using equation (3) with (T^o/2q)L.T. '^R/2q 

is indicated at the downstream end of the surface. Although the 
numerical differences between the use of- (To/ 2 q ^3 Y 'the use of 
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the sum of (^'Vo/2qjj^ rj, and '^^2q are indeed small in the case 

illustrated, the concept of diminishing skin-friction coefficient with 
increasing velocity-profile shape factor H as indicated by ('^c/2i)L,T. 
seems more logical than the concept indicated ty ^TQ/2q)g y that the 

skin-friction coefficient is independent of the velocity-profile shape. 

On this basis, the assimiption has been made that the u.se of equations (3) 
and (5) with ('i^o/2q)L.T. ''Tl/2q is more nearly correct. In 

subsequent reference to these equations by number, the use of 
(T:o/2q)L,T. '^u/^q is intended. 


Boundary-Layer Shape Parameter 

The equation given in reference 1 for the rate of growth of the 
boundary-layer shape factor is 


6 ^ = Qh.680(E-2.97^) 
dx 


_ e ^ 2q _ 2 .o3^(H - 1.286)1 

q dx TQ 


( 6 ) 


Unlike the momentum equation, this expression is basically empirical; 
moreover, it was developed for two-dimensional application and the 
correct procedure for adaptation to conical-diffuser flow is obscure. 
Attempts to use equation (6) for the calculation of diffuser flows were 
generally discouraging. An example of such a case with two-dimensional 
thickness parameters is shown in figure 3. The calculated gradients 
strongly exceed the experimental gradients and the calculated values 
of .H diverge from the experimental curve hy a steadily increasing 
margin. This example is typical of results obtained in analyses of a 
number of similar flows. In some instances, adjustments could be made 
in the numerical coefficients of equation (6) which would improve the 
calculation for a given example of diffuser flow without perceptibly 
altering the characteristics of the equation as applied to the airfoil 
flows for which it was originally developed. These cases were special, 
however, and no single set of coefficients was found that offered 
substantial improvement in some examples without detriment in others. 

Because of the difficulties and ambiguities in the use of 
equation (6), a comparable equation was sought through analysis of the 
kinetic-energy equations of the boundary layer. As a result of this 
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investigation, the following equation, for which the derivation for two- 
dimensional flow is presented in appendix A, was obtained: 


fi , ■ 1 ^ H( , 3H - 1)(H - IJ , ^ h(3H - I) - ^ - l)g (7) 

dx q.dx 2 2q ^ ' 2q 2 ^ ' 


In this equation the term Tg/2q is defined as 




T ^ , 
— ^ dy 

2 q 


( 8 ) 


Inherent in the derivation of equation (7) is the assvmption that the 
actual-boundary-layer velocity distribution may be approximated with 
acceptable accuracy by power profiles. An essentially identical equation 
in more general form and derived by a more general process has been 
presented by Tetervin and Lin (reference 12). A comparison of 
equation (7) with the comparable equation of -reference 12 is given in 
appendix B. Tetervin and Lin (reference 12) point out that the equation 
is equally applicable to three- and two-dimensional flows. 

The opportunity for determining whether equation (7) yields values 
of dH/dx in agreement with experiment when a chosen skin-friction law 
and measured values of Vs/2q are used is seriously restricted because 
of the paucity of data on the shear-stress distribution in a turbulent 
boundary layer. In reference 5, however, are found hot-wire data from 
which the values of Ts/2q may be evaluated for several flow conditions. 
When the experimental values of Tg/2q were used together with 
(T:o/2q)j^^(ji^ , equation (7) yielded shape-factor gradients in good agree- 
ment with those observed experimentally. This agreement is illustrated 
in figure U, in which the solid curve gives the shape-factor gradient 
derived from a curve faired through the experimental data, and the points 
indicated are calculated as just described. 

Equation (7) as derived and applied in the present analysis is based 
upon mean velocities and neglects tiiose- terms which would arise from 
turbulent normal stresses. The order of agreement between values from 
equation (7) and from experiment has been such as to justify this 
omission in calculating the growth of the shape factor for all conditions 
investigated. Subsequent to the preparation of the current analysis, 
other research workers at the-Langley Laboratory have supplied an 
ebctended form of the energy equation incorporating turbulent-normal- 
stress terms. Because both the momentum and kinetic-energy equations 
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are derived from the Prandtl boundary-layer equation, the inclusion of 
the Reynolds normal-stress term in the momentum equation seems to 
require its use in the kinetic-energy equation. It has been found, 
however, that, for the data of reference the net change in calculated 
rate of growth of shape factor arising frcm the inclusion of these 
additional terms is indeed quite small and in the present investigation 
they are omitted. Appendix C gives the details of this investigation. 


CORRELATION OF STRESS COEFFICIENTS 


If equations (3) and (7) are to be used for calculating boundary- 
layer development from known values of displacement thickness, momentum 
thickness, and shape factor and from stream conditions outside the 
boundary layer, it is necessary to express the three stress coefficients 
To/2q, T:s/2q, and Tjj^/2q directly or indirectly in terms of the 

assumed basic variables ^nondimensional pressure gradient ^ 

boundary-layer shape factor H, and boundaiy-layer Reynolds number Rej . 
The fundamental correctness of such a representation, however, is * 
doubtful. Dryden (reference 13) has pointed out that the turbulent 
fluctuations and turbxilent shear stress at a point in the boundary layer 
cannot be regarded as directly related to the mean speed and its deriva- 
tives at that point. The possibility exists, however^ that such 
representations, although fundamentally inexact, may be adequate for the 
purpose of calculating momentum-thickness and shape-factor gradients from 
equations (3) and (7). Furthermore, it is hoped that a more successful 
approximation for a coefficient might be obtained by an integration 
across the boundary layer than would be possible ty representing the 
value of the integrand at any local point within the boundary layer. In 
view of this possibility of adequacy for the purposes intended and the 
obvious utility of such representations if found, an investigation to 
find such representations was undertaken. 

With regard to the correlations that follow it is emphasized that 
these correlations represent those which have given the best results in 
the simplest form for the data that have been examined to date. They 
are strictly empirical and are regarded as interim representations only. 
As more data covering a broader range of conditions become available, it 
will be necessary to revise these correlations not only for greater 
numerical accuracy but also for improved correctness of representation. 

Skin-friction coefficient (T^o/2q)L.T. •“ Ludwieg and Tillmann 
(reference 3) from their experiments with the hot-plate technique 
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reported their findings in the- form of the following equation for the a 

skin-friction coefficient" in terns of the Reynolds number Hg and 
shape factor H: 


0»123 

2qyL.T. “ j^^O.268 ^q0.^78H 


( 9 ) 


This equation is plotted in figure ^ for several values of H. The 
scope of the experiments in reference 3 was insufficient to create a 
feeling of confidence that other variables, not considered in the 
analysis, may not have influenced the results. In view of the good 
results obtained in the present investigation by using this represen- 
tation, however, and the- lack of definite evidence to the contrary, the 
representation of skin-friction coefficient from reference 3 has been 
adopted for the present analysis. 


Dissipation coefficient Tg/2q.- A representation of Tg/2q as a 


function of - H, and Rg was obtained from a detailed study of 

data representing a large variety of combinations of flow conditions and 
boundary-layer parameters. This representation appeared adequate for all 


values of e The only directly measured values of ^2, avaUable 

were those of reference S. For other combinations of flow conditions 
from references 1, 8, and 9 and. data from conical-diffuser tests, 'tg/2q 
was obtained by solution of equation (7) for experimental values 


of- - — 4^ H, Rq, and 0 4^. The correlation thus obtained is pre- 
q dx^ ^ dx 

sented in figure 6. 


Over the pressure-gradient range up to - ^ = 0.008 the data from 


independent sources were sufficient"^io give reasonable confidence in the 
accuracy of the representation. The data that"vrere investigated at 
higher pressure gradients, although insufficient to establish the curves, 
were not oontradictory to the extrapolatiocrshown by the dashed lines in 
figure 6. For study of flow development'at higher pressure gradients, 
these extrapolations are suggested for use until data sufficient to 
establish their position with greater accuracy are obtained. 


The precedure by which this correlation was developed started with 
extensive plotting and cross-plotting of the coefficient T:s/2q against 

the several basic parameters - — H, and 

q dx 

search of trends. By assumption and trial of tentative algebraic 


Rg, in this order, in 






K. 
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representations of such apparent trends as could be detected, a semblance 
of order was achieved. A linear variation with the nondimensional 
pressure gradient was assumed, and an apparently systematic increase 
in Tg/2q with increasing H which appeared to be applicable over much 
of the pressure-gradient range was found. This result was cross-checked 
by giving an algebraic representation to the variation of Tg/2q with 
shape factor in order to explore the variation with nondimensional 
pressure gradient. With these general relations blocked out, the 
algebraic representations were modified and refined to incorporate the 
characteristics of low-pressure-gradient boundary-layer development with- 
out substantially altering the relationships determined at larger 
pressure gradients. In these relations, the Reynolds number factor was 
put in a form parallel <to that of Ludwieg and Tillmann (reference 3) 
for consistency. Because an orderly interrelationship of the. several 
variables developed in the process, it was possible to effect appreciable 
improvement in the representation by readjustment of the various 
constants. In graphing the final equation the Reynolds number term and 
zero-pressure-gradient - shape-fact.or term were grouped with the coef- 
ficient T:s/2q for clarity of representation. 

The order of dispersion of the data from the correlation finally 
adopted is illustrated in figure 7. In this figure the ratio of the 
experimental value of Tg/2q to the correlation-curve value is plotted 
against the nondimensional pressure gradient for five groups of data dis- 
tinguished by the ranges of shape-factor values into which the data fall. 

Turbulent-normal-stress coefficient 'TR/2q .- Values of Tjj/2q 
calculated from hot-wire data were obtained from reference ^ and 
additional values for other combinations of flow conditions and bovindary- 
layer parameters necessary for the achievement of a correlation were 
obtained from solution of equation (3) or (^) by use of data from 
references 1, 8, and 9 and frcm conical-diffuser tests. The correlation 

obtained is given in figure 8. The coordinate 0 although a derived 

quantity and not one of the basic parameters, is, however, a function of 
the basic variables. Therefore || is an implicit function of the 

basic variables. This parameter was adopted after a study of plots 
against a number of variables and combinations thereof. Although 
suggestions of trends appeared in several of the plots, llie dispersion 
of the data, wide as it is, was definitely less in the plot finally 
selected. The degree of certainty associated with the values of the 
various data points varies considerably and depends upon the method of 
determination and the test conditions. Consideration was given to this 
circumstance in constructing the curve given for the coordinates 
selected. Some of the dispersion of the points in figure 8 may arise 
from inability to fix the' value of Tjj/2q to the desired degree of 
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accuracy ty the- method employed. It is equally probable, however, that^- 
the dispersion shown may be indicative of the absence of a significant 
(jorrelating factor. 


RESULTS AND DISCUSSION 


The practical application of either equation (3) or (5) and 
equation (7) to turbulent-boundary-layer development must be regarded 
as essentially empirical. The utility of the equations is dependent on 
the success achieved in attaining correlations for the factors To/2q, 
Tg/2q, and Tj^/2q, and the degree of--success is measured by the accuracy 
obtained in applying the equations. In the figures which follow 
(figs. 9 to 19-)j "the accuracy of the analysis, in which the empirical 
representation of the-several coefficients is treated as an integral 
part, is demonstrated by comparisons between experimental- and calculated 
values of momentum thickness 6 and shape factor H for a represen- 
tative variety of examples. These demonstrations are made by two 
procedures. In the first of these procedures, the gradients of 9 
and H are calculated from equation (3) or {$) and equation (7) from 

experimentally determined values of - — H, and Ra at each boint 

q dx 

of measurement. These gradients are plotted against distance along the 
surface and integrated by means of an integraph to obtain curves of 0 
and H. These curves are intended to show the general nature of.. local 
agreement between calculated and experimental values in a manner as 
free as possible of errors arising from inaccuracies in upstream regions. 
These curves are identified by the phrase "local check." The-second 
procedure for examining the accuracy of the method is to calciilate the 
development of 9 and H as a function of longitudinal distance ty a 
step-by-step integration in which only the-stream conditions at the 
edge of -the boundary layer and the initial values of 0 and H are 
used. Curves of. 9 and H as a function of x or X, obtained by 
this- process, are identified as "predicted" curves and are intended to 
show by comparison with the "local check" curves the extent to which 
cumulative errors of integration affect -the result. 


Two-dimensional flow .- The first of these examples (fig. 9) is for 
the two-dimensional data of reference The agreement obtained is 
excellent by either method for both 0 and however, it should be 
borne in mind that these data, which were much more detailed than ary 
of the data frcan other sources, were drawn upon heavily in formulating 
the analysis. 


Figure 10 shows results as applied to the data of reference-9 
obtained from tests on a flat-plate with a boundary layer of the same 
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general size as that of reference $ but generated in a channel differing 
somewhat in geometric detail. As the end of the channel is approached, 
the calculated values of 9 deviate from the experimental values. Fig- 
ure 11 shows results as applied to the very similar tests of Wieghardt 
(reference 8 ). 

The results applied to additional data from reference 9, where the 
pressure gradient on the flat plate was adverse to only a very slight 
degree, are presented in figure 12. Agreement in H is especially 
good; the calculated values of momentum thickness 0 are slightly above 
the experimental values at the downstream end. In this example, the 
contribution of the pressure gradient to the rate of change of shape 
factor is very small, and it is to be expected that the Reynolds 
number Eq should be the dominant influence in determining the value of 
shape factor H. That this presumption appears to be valid is illus- 
strated in figure 13 in vdiich the experimental values of shape factor 
from figure 12 are plotted against the Reynolds number R 9 together With 
a curve obtained by a fairing of the experimental data taken on a flat 
plate at zero pressure gradient, both those of Ludwieg and Tillmann (ref- 
erence 3) and also those of Schultz-Grunow (reference lix). 

As a final case of two-dimensional flow, figure lij. shows the 
calculation made for the data from the example of boundary-layer develop- 
ment over a two-dimensional airfoil, cited in reference 1. Here again 
the agreement is generally good, with the notable exception of the last 
point which immediately preceded separation. In this specific example, 
the result obtained by the present method is unsatisfactory and not so 
good as that shown for this example in reference 1 . 

Conical-djf fuser flow .- Four examples of boundary-layer flow in a 
conical diffuser are shown in figures l5^ I 6 , 18, and 19. The diffuser 
for these examples was a 12° cone with a 2:1 outlet- tp-inlet area ratio 
and a 21-inch-diameter inlet. Calculations were started at the first 
measuring station downstream of the faired junction between the cylin- 
drical inlet duct and conical section. Calculations upstream of this 
point, although considered generally satisfactory, are not offered as 
evidence of the accuracy of the analysis because of the insufficiency of 
experimental data in this region where the pressure changes most rapidly. 
The first of this group of figures (fig. 1^) is for an inlet boundary 
layer having a momentum thickness of 0.00217 foot and a shape factor 
of 1.226 at a mean Mach number in the approach duct of 0.27. The second 
(fig. 16 ) is for the same Mach number with an inlet-boundary-layer 
momentum thickness of 0.0112 foot and a shape factor of 1.220. For both 
of these figures the agreement between t he oiy and experiment is quite 
good except at the final points. The station at viiich the final points 
were obtained actually was in the upstream end of the cylindrical 
discharge duct rather than in the cttf fuser proper. Because of this' 
circumstance, the validity of using these points to assess the accuracy 
of the procedure is questionable. 


-1 
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The conical-diffuser flow illustrated in figure 16 is that used in 
figure 3 to illustrate the inadequacy of equation (6). The results 
obtained with equation (7) plotted to the same coordinates as those of 
figure 3 are shown in figure 17. The improvement in agreement as 
obtained in equation (7) is evident from a comparison of figures 3 and 17. 

Although the present analysis is for incompressible flowj the 
results appear to be acceptable- for moderately high speeds, as illus- 
trated by figure 18, in which calculations are presented for an initial 
boundary-layer momentum thickness of 0.001i;8 foot, a snape factor 
of 1.225, and a mean inlet Mach number of- 0.^3. The principal discrep- 
ancy is at the last— points. In figure 19, calculations are shown for an 
initial momentum thickness of 0.0112 foot, a sh^pe factor of 1.223, 
and a mean inlet Mach number of 0.5U. 


Separation .— It is not the purpose of this paper to propose new 

criteria for the calculation of the point of separation of the boundary 
layer nor to consider the relative merits of such criteria as may have 
been proposed by others. It is of interest to note, however, that in the 
process of making the calculation of the boundary-layer development ‘for 
the data of reference 5 the trend of the calculated skin-friction- 
coefficient curve was observed to be troward zero at a point at which 
separation had been established by direct measurements in the originsil 
experiment. This trend is illustrated in figure 20. Although the 
Ludwieg and Tillmann correlation (reference 3) makes the skin-friction 
coefficient L.T. asymptotic to zero as the shape factor H 

becomes infinite, at the last point at which measurements were made 
prior to the qccurrence of separated flow, it appears that the rate of 
increase of the shape factor H with distance is itself increasing at 
a rate sufficient to give the curve of skin-friction coefficient 
^To/2q)]^_,j.^ a trend suggestive of a finite intercept on the X-axis. 


APPLICATION 


The preceding discussion has been concerned with the accuracy -of the 
analysis which was examined by application to a representative groi^) of 
examples by two methods. Inasmuch as the principal anticipated utility 
of the current analysis would be the prediction of. the boundary-layer 
development; the accuracy of the analysis, as. demonstrated by the second 
procedure (perf onning a stepwise integration using only the initial 
values of 6 and H and stream conditions outside the boundary layer), 
is of particular interest. The necessity of repeated solution of 
either equation (3) or (5) and equation (7), however, renders this step- 
wise calculation especially inconvenient. In order to facilitate 
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execution of this procedure j the terns of either equation (3) or (3) 
and equation (7) have been grouped and graphed in a convenient manner. 

The analysis as so organized is given by the following equations together 
with associated figures: 

e M = I + II (10) 

dx 

For two-dimensional flow: 


dB 

dx 


= in + IV + V 


( 11 ) 


For conical diffuser flow; 


^=III + IV+V--^^ 
dx D dx 


( 12 ) 


Algebraic expressions corresponding to the graphs of terms I, II, III, 
IV, and V are given in appendix D. 

The terms I and n in equation (10) are obtained from figures 21 
and 22, respectively, as functions of the nondimensional pressure 

gradient ~ q ^ the shape factor H. Term I is read directly from 
figure 21. Term II is obtained by dividing the ordinate of figure 22 
by Although the curves of figure 22(a) appear to have 

e dq 

a common intercept at q ^ scaae displacement of the intercept 

actually occurs for different values of H, and this displacement is of 
importance vdien developments at very small pressure gradients are to be 
calculated. An enlargement of the graph for small pressure gradients 
is given in figure 22(b), and the zero-pressure-gradient intercepts are 
plotted as a function of shape parameter in figure 22(c). Term III for 
equations (11) and (12) is given in figure 23 as a function of the 

nondimensional pressure gradient “ 'I shape factor H. Term IV 

is ('’^o/2q^L^ip^ and is represented in figure 3 as a function of the 

Reynolds number Rq and shape factor H. Term V is TR/2q and is 
represented in figure 8 as a function of the shape-factor-growth coef- 
ficient e ^ which is made available by the solution of equation (10). 

For calciilating the boundary-layer growth in a conical diffuser the duct 
dimensions must be known and equation (12) should be Used. 
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CONCLUSIONS 


A procedure based on IJrie kinetic-energy equation and an extended' 
form of the momentum equation has been devised for calculating the 
development of turbulent boundary layers, in adverse pressure gradients, 
for that class of flows for which the fluid density at all points and 
total pressure outside the boundary layer are invariant. In the develop- 
ment of this procedure an effort was made not only to arrive at an ana- 
lytical form that would allow examination of the significance’ of the 
physical quantities involved but^^so to achieve a high degree of con- 
sistency with the more recent -results of turbulent-boundary-layer 
research. 

Correlations, which are essential to the execution of the method, 
of several elemen'ts of the basic equations in terms of quantities 
regarded as controlling parameters are presented. These correlations 
are strictly interim empirical relationships both as to numerical values 
and as to the variables involved. It is to be expected that, with the 
acquisition of more data covering a broader range of conditions, improve- 
ments will be made— in both the functional nature and the accuracy of the 
correlations. A comparable and associated refinement of the basic 
equations is likewise anticipated as more data are obtained. 

Predictions, by this method, of turbulent-boundary-layer development 
are compared with experimental results from several sources .for a number 
of cases of flow over flat plates and airfoils and in conical diffusers. 
In the range of boundary-layer flow short cf separation, the agreement 
of calculated values-with experimental values was, in most cases, quite 
satisfactory; in some instances, however, definite~disagreement between 
the calculated and experimental results was noted. It is believed, 
however, that good agreement has been obtained in enough instances to 
Justify conttouation of effort along the present lines, particularly 
with respect to improvement of the correlations and refinement -of the 
equations . - 


Langley Aeronautical Laboratory 

National Advisoiy Committee for Aeronautics 
Langley Field, Va., April 26, 19^1 
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APPENDIX A 

DERIVATION OF EQUATION FOR RATE OF GROWTH OF BODNDARI-LAYER 
SHAPE FACTOR (EQUATION (7)) 


The two-dimensional equation for the loss of kinetic energy per 
unit length, normal to the direction of flow, is given in reference 7 as 




- I dy 


(Al) 


The kinetic-energy thickness of the boundary layer is defined as 



If equation (Al) is divided by U and the definition for the 
kinetic-energy thickness 63 is substituted in the right-hand side of 
the equation, the following equation results; 



^ dy = — 
3y U? 


du3d3 

dx 


q 


By substitution of the definition for. Ts/2q 


(A2) 


P ^ 

2q Jo 2q By 


equation (A2) becomes 
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Sy differentiating the right-hand side of equation (A3) 


Is = tt3 ^^3 . 3 W 

2q 2\]3 ^ ^ dx 


and substituting the relation 


pU 


dU ^ ^ 
dx dx 


the following equation is obtained; 


Is = + 1 £2 

2q 2\dx 2 q dx/ 


(Ah) 


If the shape of the velocity profile is assumed to be represented 
by the relation 


!-( 5 ) 


1/n 


where 


n = 


H - 1 


it can be shown that kinetic-energy thickness 63 may be represented 

by 


(A^) 


c- _ n ha 

^3 = e 3r=rT 
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By transposing terms and substituting the representation given 
for 63 , the following equation is obtained from equation (AU): 


^53 ± L _Uh\= 9 _ 1 ^ 

dx * dx ^ 3H - ly 2q 2 q dx 

Differentiating the left-hand side of tliis equation and transposing 
terms yields an equation of the following form: 


e ^ = H(3H - 1) 


dB 

dx 


3H - i f^ '^s _ 3 63 
Hh ^ ~ 2 q dxy 


(A6) 


Substituting the momentum equation (equation (1)) and the definition for 
63 (equation (A^)) into this relation yields the kinetic-energy equation 
for the rate of growth of boundary-layer shape factor H in its final 
form for the present analysis 
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APPENDIX B 

COMPARISON OF EQUATION (7) OF THE PRESENT ANALYSIS 
WITH EQUATION (13) OF REFERENCE 12 


Equation (13) in reference 12 is given as 


0 - I ^ H(H - 1)(3H - 1) + (3H - r)(H + 


3H - 1 




2 Jo ^ 


pU2 


12 (H + 1)(3H - 1) 
U T 


(Bl) 


with the following accompanying definitions; 




Vq is the velocity perpendicular to the Wall at the wall (y ■ 0), 
and p is an exponent related to th,e profile shape. Substitution 
of these definitions in equation (Bl) and use of the relation 

U dx ” 2q dx 
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gives the following equation: 


dH _ 0 dq H(3H - 1)(H - 1) , . 

dx ■ q dS 2 ^ ^ ^ 


-I 

'^0 (3H - 1)2 p-*- u '^O ^ y , "^0 (H + 1)(3H - 1) 

5q 2 Jo D^<»6*T -E 

6 

The relation 



(B2) 


may be nondimensionalized in the following manner: 


r H d I = 
Jo s 
6 


pi At 

T U/U 

Jo •'0^ 



Substitution of this relation in equation (B2) yields an equation of the 
following form: 


e £ . - i ^ h(3k-i)(h-i) , ^ „ _ _ 

dx q dx 2 2q ^ 

'^0 (3H - 1)2 T ^ , y , ^0 (H + 1)(3H - 1) 

2q 2 Jo E 


(B3) 


Substituting the definition for Tg/2q 
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:mto equation (B3) and noting that the term Vq/U becomes zero for the 
solid-wall case being considered herein yields 

9 dq H(3H - 1)(H - 1) , .n ^s (3H - l)^ 

dx " q dx 2 2q ^ 2q 2 

The result obtained is the same as equation (?) of the present analysis. 
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APPENDIX C 

EXTENDED ENERGY EQUATION 


The energy equation as extended to include turbulent-nomial-stress- 
coefficient terms is: 


0 ^ i dq H(H - 1)(3H - 1) ^ 2° Hf3H - 1) - — 

dx q dx 2 2q 2q 2 


^ hoh - 1) - (01) 


in which, as before, 


Is = Jl P a(pu 

2q"2qjQ dx 


and the new term Tjj*/2q is defined as 


V ^ ^ u a.(,pu', ^) ^ 

2q 2qJo U dx ^ 

It may be noted that equation (Cl) differs from equation (7) by the 
addition of the terms in and Tjj»/2q. The sum of these two 

terms, divided by 0, therefore, constitutes the correction to the calcu- 
lated shape-factor gradient which is to be made to provide for turbulent 
normal stresses. Figure 2k shows the experimental data of reference 5 
and the shape-factor development calculated from equation (7), trans- 
ferred from figure $• (The curve of dH/dx, which was integrated to 
obtain the curve of H in figure 9, has been added for more complete 
comparison.) The dashed curve of figure 2l;(a) has been obtained by 
adding the turbulent-normal-stress corrections to the solid curve, and 
the integration of this dashed curve constitutes the dashed curve of H 
against x in figure 2i;(b). In general, the departures, either in rate 
of change or value of the shape factor, are not great. In interpreting 
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figure 2k, the significant item is the deviation of the dashed curves 
from the solid curves, not the general position of either relative to 
the data, because the empirical representation of Ts/2q was established 
without regard to the turbulent-normal-stress terms. 
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ALGEBRAIC EXPRESSIONS FOE TERMS OF EQUATIONS (10), (11), AND (12) 


The algebraic expressions corresponding to the graphs of terms I, II, m, IV, and V are 
as follows; 


Term I 


e dq (H - 1)(3H - 1) 
q dx 2 


Term II x Rg 


0.268 - _ i 
2 


X 10-0.678H |o.02U6(H - 1) - 


I ^[i67.2(H - 1.2?)0*^35 _ 35.1^] [h(3h ^ 0.9) 



Term III 



Term IV = 


0.123 

0.268 jq0.676h 


/ dTT 

Term V = 0.0011(9 — + 0.1876 — 

\ dx/ dx 

It should be particularly noted that these algebraic representations constitute a mathe- 
matical fitting of the correlations obtained and shotild not under any ciiTCtimstances be regarded 
as having any other significance or applicability outside the range of the curves shown. 


ro 
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(b) Momentum tblckness. 

Figure 1.- Comparison of calculated and experimental momentum- thickness 

gradient and momentum thickness. 
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Figure 2.- Variation with, longitudinal distance of skin- friction 
coefficient and turhulent-norml-stress coefficient. 
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(a) Shape- factor gradient. 
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(h) Shape factor. 

Figure 3.- Congparison of calculated (by method of reference 1) and 
eiqperlmental shape- factor gradient and shape ftictor for conical- 
diffuser flow for Mach number 0.27. 
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PiguiB 0 .- Correlatioa of turbiilent-nomial-stress coefficient 
repreBented as a function of the nondlmenalonal shape-factor 

gradient 0 
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(a) Local check. 

Figure 9.- Variation with longitudinal distance of eaqperlmental and 
calculated values of momentum thickness and shape factor. 













(a) Local check. 


Figure 10.- Variation with longitudinal distance of experimental and 
calculated values of momentum thickness and shape factor. 
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("b) Predicted. 
Figure 10.- Concluded. 





(a) Local check. 


Figure 11.- Variation with longitudinal distance of experimental and 
calculated values of momentum thickness and shape factor. 
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^■b) Predicted. 


Figure 11.- Concluded. 




(a) Local cLeck. 

Figure 12.- Variation with longitudinal distance of experimental and 
calculated values of momentum thickness and shape factor. 
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(a) Local check. 


ik. ~ Variation wi'fch loagi'fcad.iiial d.ls'taace of ©s^erimea'fcal and. 
calculated values of momentum thickness and shape factor. 








(b) Predicted. 
Figure Concluded 
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(■b) Predicted. 
Figure l6.- Concluded. 
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(a) - 1 1| = 0 to O.OlOj H = 1.286 to 2.00. 

S^gure 22.- Term II x EgO .268 

of equation (lO) expreaaed as a function 
of the nondinensional pressure gradient - ^ ^ for constant values 
of shape fhctor H. 
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Wgure 23 .- term IH of eiuotlone (U) and (12) expi^ased aa a function 

of the nondi.KnBlonal preaauxo gradient -2^ for conatant talnee 

of ahape f&.ctor H. 












